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Neutrino oscillations can explain the observed motion of pulsars. We show that two different models 
of neutrino emission from a cooling neutron star are in good quantitative agreement and predict the 
same order of magnitude for the pulsar kick velocity, consistent with the data. 
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I. INTRODUCTION 



We recently suggested jjJH that the observed proper 
motions of pulsars || may be the result of neutrino os- 
cillations in the hot neutron star born in a supernova ex- 
plosion. Neutrinos are the only significant cooling agents 
during the first 10 seconds after the onset of the super- 
nova, and they carry away most of the energy liberated 
in the gravitational collapse, ~ 10 53 erg. A 1% asymme- 
try in the distribution of neutrinos can account for the 
measured pulsar velocities ~ 500 km/s. 

Neutrino oscillations in medium are affected by the po- 
larization effects due to magnetic field Q . Consequently, 
the strong magnetic field of a neutron star can affect the 
depth at which the neutrino conversions take place. The 
points of the Mikheev-Smirnov-Wolfenstein resonance lie 
on a surface that is not concentric with the star. For a 
certain range of neutrino masses, this creates an asym- 
metry in the distribution of momenta of the outgoing 
neutrinos. The <-» v e oscillations as well as the 
sterile-to- active neutrino conversions [Q] , can give the pul- 
sar a kick consistent with the observation provided that 
the magnetic field inside the star is of order 10 14 — 10 15 G. 
Other types of neutrino conversions could also produce a 
similar effect ||. 

Our explanation |l]J^| uses the fact that different neu- 
trino species have different opacities in nuclear matter. 
The charged-current interactions are responsible for the 
difference in the mean free paths of v e and v^^. They 
are also the cause of the matter-enhanced v e <-» v^ r con- 
versions. In the case of sterile neutrinos, both charged 
and neutral currents contribute to the difference in the 
opacities. The analyses of Refs. QJ^H was based on a 
simplified model which, as we will see, gave the right 
order-of-magnitude estimate for the pulsar velocity. 

Recently, a different model for the neutrino emission 
was used to calculate the kick HI due to the active neu- 



trino oscillations]]. The result obtained in Ref. Iq] was in- 
correct because the neutrino absorption v e n — > e~p + was 
neglected and also because the different neutrino opaci- 
ties were assumed to be equal to each other. We empha- 
size that in the absence of charged-current interactions 
the kick from the active neutrino oscillations Q should 
vanish. 

We will show that, after the charged-current interac- 
tions are included, the two models are, in fact, in good 
agreement, as they should be. 

II. HARD NEUTRINOSPHERES 

Our calculations employed a model with a sharp 
neutrinosphere, such that the neutrinos of a given type 
were assumed trapped inside and free-streaming outside. 
Several subsequent analyses || used the same model to 
calculate the kick to the pulsar from the neutrinos whose 
opacities changed on their passage through matter due 
to oscillations. 

In this model one assumes that the luminosity of the 
emitted neutrinos obeys the Stefan-Boltzmann law and 
is proportional to T A . 

This model is admittedly simplistic. However, it ap- 
pears to work well insofar as predicting the order of mag- 
nitude of the kick. 



III. SOFT NEUTRINOSPHERES 

Let us now consider the Eddington model for the at- 
mosphere which was used by Schinder and Shapiro to 



J The active-to-sterile neutrino conversions of Ref. || are as- 
sumed to take place at much higher densities than those dis- 
cussed in Ref. M. 
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describe the emission of a single neutrino species. We 
will generalize it to include several types of neutrinos. 

In the diffusion approximation, the distribution func- 
tions / are taken in the form M: 



eg 



A, 9m 



(1) 



where / e<? is the distribution function in equilibrium, Aj 
denote the respective opacities, m is the column mass 
density, m = J pdx, £ = cosa, and a is the normal angle 
of the neutrino velocity to the surface. At the surface, 
one imposes the same boundary condition for all the dis- 
tribution functions, namely 



/„ 4 (ro,O=0, for£<0, 
f Vi {m,t) = V eq , for£>0. 



(2) 



However, the differences in Aj produce unequal distribu- 
tions for different neutrino types. 

Generalizing the discussion of Refs. [BQ to include six 
flavors, three neutrinos and three antineutrinos, one can 
write the energy flux as 



F = 2n E 3 dE / £d£ Y7/„ 4 + 



(3) 



We will assume that A, = a!?\e 2 /E 2 ). 

We use the expressions for from equation ([!]). 
Changing the order of differentiation with respect to m 
and integration over E and £, and using the fact that 
f eq is isotropic, we arrive at the result similar to that of 
Ref. f): 



F 
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dT 2 
dm 



(4) 



The basic assumption of the model is that flux F 
is conserved. In other words, the neutrino absorptions 
v e n — > e~p + are neglected. Since the sum in brackets, 



as well as the flux F are treated 
respect to m, one can solve for T 2 : 



as constants with 
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Swapping the two flavors in equation (||) leaves the 
temperature unchanged in the Eddington approximation. 
Hence, neutrino oscillations do not alter the temperature 
profile in this approximation^. 



2 It was also assumed in our earlier calculations that 
neutrino oscillations do not have a significant effect on the 
temperature profile. We then took the spectral temperature 
of free-streaming neutrinos to be equal to matter temperature 
at their respective neutrinospheres. 



We will now include the absorptions of neutrinos. 

Some of the electron neutrinos are absorbed on their 
passage through the atmosphere thanks to the charged- 
current process 



e p 



(6) 



The cross section for this reaction is a — 1.8G 2 F E 2 , where 
E v is the neutrino energy. The total momentum trans- 
fered to the neutron star by the passing neutrinos de- 
pends on the energy. 

Both numerical and analytical calculations show that 
the muon and tau neutrinos leaving the core have much 
higher mean energies than the electron neutrinos pH. 
Below the point of MSW resonance the electron neutrinos 
have the mean energies « 10 MeV, while the muon and 
tau neutrinos have energies « 25 MeV. 

The origin of the kick in this description is that the 
neutrinos spend more time as energetic electron neutri- 
nos on one side of the star than on the other side, hence 
creating the asymmetry. Although the temperature pro- 
file remains unchanged in Eddington approximation, the 
unequal numbers of neutrino absorptions push the star, 
so that the total momentum is conserved. 

Below the resonance E Ve < E Vt . Above the reso- 
nance, this relation is inverted. The energy deposition 
into the nuclear matter depends on the distance the elec- 
tron neutrino has traveled with a higher energy. This 
distance is affected by the direction of the magnetic field 
relative to the neutrino momentum. 

We assume that the resonant conversion v e <-> v T takes 
place at the point r = ro + 6 ((f)); 5 ((f)) = Sq cos cf>. The 
position of the resonance depends on the magnetic field 
B inside the star [BJ: 



So 



2ir 2 



'dN e 
dr 



(7) 



where N e — Y e N n is the electron density and /x e is the 
electron chemical potential. 

Below the resonance the r neutrinos are more ener- 
getic than the electron neutrinos. The oscillations ex- 
change the neutrino flavors, so that above the resonance 
the electron neutrinos are more energetic than the r neu- 
trinos. The number of neutrino absorptions in the layer 
of thickness 26(<f>) around r depends on the angle <f> be- 
tween the neutrino momentum and the direction of the 
magnetic field. Each occurrence of the neutrino absorp- 
tion transfers the momentum _E„ c to the star. The dif- 
ference in the numbers of collisions per electron neutrino 
between the directions 4> and — is 

Ak e /E Uc = 2 5{cp) N n [o-(Si) - o-CEa)] (8) 

= 1.8 G 2 F [E 2 -E 2 ]^^-h Nc cos^, (9) 

where h Ne = [d(\n N e ) / dr]^ 1 . 

We use Y e w 0.1, E x w 25 MeV, E 2 ~ 10 MeV, 
fi e ps 50 MeV, and /iat c k 6 km. After integrating over 
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angles and taking into account that only one neutrino 
species undergoes the conversion, we obtain the final re- 
sult for the asymmetry in the momentum deposited by 
the neutrinos: 



Afc 



= 0.01 



B 



2 x 10 14 G' 



(10) 



which agrees with the earlier estimates^] EIH. 

Neutrinos also lose energy by scattering off the elec- 
trons. Since the electrons are degenerate, the final-state 
electron must have energy greater than /i e . Therefore, 
electron neutrinos lose from 0.2 to 0.5 of their energy per 
collision in the neutrino-electron scattering. However, 
since N e <C N n , this process can be neglected. 

One may worry whether the asymmetric absorption 
can produce some back-reaction and change the temper- 
ature distribution inside the star altering our result (|To|). 
If such effect exists, it is beyond the scope of Eddington 
approximation, as is clear from equation (|^). The only 
effect of the asymmetric absorption is to make the star 
itself move, in accordance with the momentum conserva- 
tion. This is the origin of the kick ([To|). 

Of course, in reality the back-reaction is not exactly 
zero. The most serious drawback of Eddington model, 
pointed out in Ref. 0, is that diffusion approximation 
breaks down in the region of interest, where the neu- 
trinos are weakly interacting. Another problem has to 
do with inclusion of neutrino absorptions and neutrino 
oscillations M. However, to the extent we believe this 



approximation, the pulsar kick is given by equation (10) 



IV. COMPARISON OF THE TWO MODELS 

We have juxtaposed two models for neutrino trans- 
fer. One of them assumes that neutrinos are trapped 
inside a sharply defined neutrinosphere, but that they 
are free-streaming outside. In this model, the absorption 
of neutrinos outside the neutrinosphere is neglected and 
the asymmetry arises from the change in the position of 
the neutrinosphere due to neutrino oscillations. 

The other model considers the propagation of the neu- 
trinos through the atmosphere interpolating between the 
region where the neutrinos are trapped and that in which 
they are free-streaming. Here the neutron star receives a 
kick from the unequal numbers of reactions v e n — > p + e~ 
on either side. The higher-energy neutrinos spend more 
time as electron neutrinos on one side of the star than 
on the other side. Therefore, the passage of neutrinos 
creates an uneven drag to the star. 



3 We note in passing that we estimated the kick in Refs. 
assuming [i e ~ const. A different approximation, Y e ~ const, 
gives a somewhat higher prediction for the magnitude of the 
magnetic field |lo| . 



Incidentally, the Stefan-Boltzmann relation between 
the luminosity and temperature is, of course, present in 
both models. In the diffusion approximation, equation 
(|) implies that F tx T 4 at m = 0. 

In both models the kick is a manifestation of the un- 
equal neutrino opacities which are caused by the charged- 
current interactions. 

The two models are in good quantitative agreement. 



V. CONCLUSION 

The neutrino oscillations can be the explanation of the 
pulsar motions. Although many alternatives have been 
examined, all of them fail to explain the large magnitudes 
of the pulsar velocities. Parity violation effects recently 
claimed to explain the kicks |Ll|] turned out to be irrele- 
vant because they wash out in statistical equilibrium and 
produce no appreciable asymmetry in the neutrino lumi- 
nosity p"2| . Other proposed mechanisms required unusu- 
ally high neutrino magnetic moments fl3f| or some exotic 
interactions. The oscillations appear to be the only vi- 
able explanation at present. 

If the pulsar kick velocities are due to v e <-> u^ T con- 
versions, one of the neutrinos must have mass ~ 100 eV 
(assuming small mixing) and must decay on the cosmo- 
logical time scales not to overdose the Universe |l| . This 
has profound implications for particle physics hinting at 
the existence of Majorons [JL4J or other physics beyond 
the Standard Model that can facilitate the neutrino de- 
cay. 

If the active-to-sterile neutrino oscillations 0] are re- 
sponsible for pulsar velocities, the prediction for the ster- 
ile neutrino to have a mass of several keV is not in contra- 
diction with any of the present bounds. In fact, the ~keV 
mass sterile neutrino has been proposed as a dark-matter 
candidate ]l5| . 

To summarize, two different models of neutrino emis- 
sion predict the same order of magnitude for the pulsar 
birth velocity caused by neutrino oscillations in the cool- 
ing neutron star. 
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